Previous experience with ray-tracing sun-disk images on off-axis concave mirrors, led us to investigate trajectories within elliptical billiards.
In 2011, we produced the following animation, which uncovered that certain trajectories were perioic orbits, a gallery of which can be found here.
Shortly after, we created the following video to study the geometry of triangular orbits in the ellipse, where we showed the locus of the orbits’ incenters were apparently elliptical.
Then the following fortuitous sequence of events took place (corrections welcome 😄):
To summarize: indeed, the locus of triangular orbits’ incenters is elliptical as are that of many other triangular orbit centers!
Below we reproduce the original animation this time using the R ecosystem. Without loss of generality, an ellipse is shown with semi-axes of lengths a = 2, b = 1 centered at the origin (the “x”). The ellipse’s foci, at ±sqrt(a²-b²) = ±sqrt(3), are shown as black dots. A base vertex P(t) is chosen on the ellipse at [a*cos(t),b*sin(t)], for some t in [0,2π] (red dot). Via numeric optimization we compute the other two vertices, shown blue, under the constraint they be bisected by their local normals, shown as inward-pointing arrows. For a given t the incenter (shown as a green dot) is computed by intersecting any two normals. Its locus (shown as dashed green ellipse) is obtained by chaining all points under some angular step.
Notice that the two foci are always contained within the orbits, owing to the fact that via Poncelet’s theorem, the orbits’ caustics are confocal with the original ellipse [insert ref].
We also verify a fascinating property about triangular orbits in the elliptic billiard: the perimeter of all orbits is constant (though not the area). The graph below shows the perimeter and area numerically calculated for triangular orbits for all t.
On perimeter constancy, via email Prof. Tabachnikov points out that:
the perimeter of a periodic orbit is constant (for any period, starting with three) since these n-gons are extrema of the perimeter function on inscribed n-gons. This function is constant on the curve consisting of its critical points. A somewhat related fact: the centers of mass (with two possible meanings) of a Poncelet polygon move along conics, as proven in our hoax paper with Rich Schwartz.
And Prof. Koiller adds that:
constancy of length comes from Hamiltonian Dynamics. The family of triangles forms an abstract 2-torus in phase space, all trajectories having same velocity and period, hence same length.
We developed a geometric method which yields closed-form expressions for the extrema in x and y of the incenter locus. Since we know the latter is elliptical, these correspond to the length of the two semi-axes.
Namely, the method comprises the following steps:
Simulation code (in R) can be found here
© Dan S. Reznik, 2019